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ABSTRACT 



Aims. The aim of the present work is to constrain the Coma cluster magnetic field strength, its radial profile and power spectrum by 
comparing Faraday Rotation Measure (RM) images with numerical simulations of the magnetic field. 

Methods. We have analyzed polarization data for seven radio sources in the Coma cluster field observed with the Very Large Array 
at 3.6, 6 and 20 cm, and derived Faraday Rotation Measures with kiloparsec scale resolution. Random three dimensional magnetic 
field models have been simulated for various values of the central intensity Bo an d radial power-law slope n, where rj indicates how 
the field scales with respect to the gas density profile. 

Results. We derive the central magnetic field strength , and radial profile values that best reproduce the RM observations. We find that 
the magnetic field power spectrum is well represented by a Kolmogorov power spectrum with minimum scale ~ 2 kpc and maximum 
scale ~ 34 kpc. The central magnetic field strength and radial slope are constrained to be in the range (Bo =3.9 fiG; rj =0.4) and 
(B =5.4 yuG; rj =0.7) within la. The best agreement between observations and simulations is achieved for B = 4JfiG;n =0.5. 
Values of Bo >7 /iG and < 3 yuG as well as rj < 0.2 and rj > 1.0 are incompatible with RM data at 99% confidence level. 

Key words. Cluster of galaxies - Magnetic field - Polarization - Faraday Rotation Measures - A1656 Coma 



c3 1. Introduction 



It is now well established that the intracluster medium (ICM) 
of galaxy clusters is not only composed of thermal gas emitting 
in the X-ray energy band, but also of magnetic fields permeat- 
ing the entire cluster volume (see Ferrari et al., 2008 for a re- 
cent review). This is directly demonstrated by the detection of 
large, diffuse synchrotron radio sources such as radio halos and 
radio relics, in an increasing number of galaxy clusters (see e. g. 
Venturi et al. 2008, Giovannini et al. 2009). In these clusters it 
is possible to estimate the average ICM magnetic field under the 
minimum energy hypothesis (which is very close to equiparti- 
tion conditions) or by studying the Inverse Compton hard X-ray 
emission (e.g. Fusco Femiano et al. 2004). 
The ICM magnetic field is also revealed by the analysis of po- 
larized emission of radio sources located at different projected 
distances with respect to the cluster center. The interaction of 
the ICM, a magneto-ionic medium, with the linearly polarized 
synchrotron emission results in a rotation of the wave polariza- 
tion plane (Faraday Rotation), so that the observed polarization 
angle, ^obs at a wavelength A differs from the intrinsic one, int 
according to: 



¥ obs (A) = + A 2 x RM, 



(1) 



where RM is the Faraday Rotation Measure. This is related to the 
magnetic field component along the line-of-sight (B//) weighted 
by the thermal gas density (n e ) according to: 



RMoc ( n e (l)B 1 i(l)dl. 

Jlos 



(2) 
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Therefore, once the thermal gas density distribution is inferred 
from X-ray observations, RM studies give an additional set of 
information about the cluster magnetic field. This is the only 
way, so far, to study the intracluster magnetic field in clusters 
where diffuse radio emission sources are not directly observed. 
The Coma cluster magnetic field has been studied in the past 
using all three approaches mentioned above. The first investiga- 
tion of the magnetic field was performed by Kim et al. (1990). 
They analyzed 18 bright radio-sources in the Coma cluster 
region, obtaining RM maps at ~ 20" (~ 9.2 kpc) resolution 
and found a significant enhancement of the RM in the inner 
parts of the cluster. Assuming a simple model for the magnetic 
field reversal length, they derived a field strength of ~ 2 muG. 
A complimentary study was performed by Feretti et al. (1995) 
studying the polarization properties of the extended radio galaxy 
NGC 4869. From the average value of RM and its dispersion 
across the source, they deduced a magnetic field of ~6 fiG 
tangled on scales of ~ 1 kpc, in addition to a weaker magnetic 
field component of ~0.2 /uG , uniform on a cluster core radius 
scale. 

From the Coma radio halo, assuming equipartition, a magnetic 
field estimate of ~ 0.7 -1.9 yuG is derived (Thierbach et al. 
2003), while from the Inverse Compton hard X-ray emission 
a value of ~0.2 yuG has been derived by Fusco Femiano et 
al. (2004), although new hard X-ray observations performed 
with a new generation of satellites did not find such evidence 
of non-thermal emission (Wik et al. 2009 using XMM and 
Suzaku data, Lutovinov et al. 2008 using ROSAT, RXTE and 
INTEGRAL data, Ajello et al. 2009 using XMM-Newton, 
Swift/XRT, Chandra and BAT data; see Sec.[8). 
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However, the discrepancy between these values is not surpris- 
ing: equipartition and IC estimates, in fact, rely on several 
assumptions, and are cluster volume averaged estimates, while 
the RM is sensitive to the local structures of both the thermal 
plasma and the cluster magnetic field component that is parallel 
to the line of sight. Furthermore, the equipartition estimate 
should be used with caution, given the number of underlying 
assumptions. For example, it depends on the poorly known 
particle energy distribution, and in particular on the low energy 
cut-off of the emitting electrons (see e. g. Beck & Krause 2005). 
These different estimates provide direct evidence of the complex 
structure of the cluster magnetic field and show that mag- 
netic field models where both small and large scale structure 
coexist must be considered. The intracluster magnetic field 
power spectrum has been investigated in several works (Vogt 
& Ensslin 2003, 2005; Murgia et al., 2004; Govoni et al. 
2006, Guidetti et al. 2008, Laing et al. 2008). In addition, a 
radial decline of the magnetic field strength is expected from 
magneto-hydrodynamical simulations performed with different 
codes (Dolag et al. 1999, 2005; Brueggen et al. 2005, Dubois et 
al. 2008, Collins et al. 2009, Dolag & Stasyszyn 2008, Donnert 
et al. 2009a), as a result of the compression of thermal plasma 
during the cluster gravitational collapse. 

In this paper we present a new study of the Coma cluster 
magnetic field. We analyzed new Very Large Array (VLA) 
polarimetric observations of seven sources in the Coma cluster 
field, and used the FARADAY code developed by Murgia et al. 
(2004) to derive the magnetic field model that best represents 
our data through numerical simulations. 

The Coma cluster is an important target for a detailed study of 
cluster magnetic fields. It is a nearby cluster (z=0.023), it hosts 
large scale radio emission (radio halo, radio relic, bridge) and 
a wealth of data are available at different energy bands, from 
radio to hard X-rays. 

The paper is organized as follows: in Sec.|2]the properties of the 
X-ray emitting gas are described, in Sec. [3] radio observations 
are presented, and the radio properties of the sources are 
analyzed. RM analysis is presented in Sec. [4] while in Sec|5] 
the adopted magnetic field model is described. The method we 
adopted to compare simulations and observations is reported in 
Sec. [6] numerical simulations are presented in SecQand results 
are discussed in Sec. [8] Finally, conclusions are reported in Sec. 
i 

A ACDM cosmological model is assumed throughout the paper, 
with H = 71 km/s/Mpc, £2 M =0.27, C2 A =0.73. This means that 
1 arcsec corresponds to 0.460 kpc at z=0.023. 



2. Thermal component from X-ray observations 

The study of the magnetic field through the Faraday RM requires 
knowledge of the properties of the thermal gas (see Eq. [2j. This 
information can be derived from X-ray observations. In Fig. Q] 
the X-ray emission of the Coma cluster is shown in colors. X-ray 
observations in the energy band 0. 1 -2.4 ke V have been retrieved 
from the ROSAT All Sky Survey data archive. After background 
subtraction the image has been divided by the exposure map and 
smoothed with a Gaussian of <x = 60". 

The radio contours of the NVSS (NRAO VLA Sky Survey) at 1 .4 
GHz are overlaid onto the X-ray emission in Fig.Q] The location 
of the sources studied in this paper are marked with their names. 
Note that the extended radio emission of the radio halo, relic and 
bridge are completely resolved out in the NVSS image due to the 
lack of extremely short baselines. 



The X-ray emission is from thermal bremsstrahlung, and can be 
used to trace the thermal particle distribution in the ICM. The 
distribution of the gas is well reproduced by the so-called "/J- 
model" (Cavaliere & Fusco Femiano, 1976), reported in Eq.[3] 

n e {f) = n \l + , (3) 

where r indicates the radial distance from the cluster center, r c is 
the cluster core radius, no is the cluster central gas density. The 
cluster center is RA = 12h 59m 41.5s; DEC = 27° 56' 20". We 
use the values derived by Briel et al. (1992), corrected for the 
cosmology adopted in this paper. They found : 

- =0.75+0.03; 

- r c =291 + 17 kpc; 

- « =3.44+0.04 10~ 3 cm- 3 . 

3. Radio observations and images 

3.1. VLA observations and data reduction 

We selected from NVSS a sample of sources having a peak 
flux density larger than 45 mJy, located in a radius of 1° 
from the cluster X-ray center center (^ 5r c ), and which have 
indication of polarization from Kim et al. (1990). A further 
selection was performed on the basis of the position of the 
sources. Observations have been performed at the VLA on 7 
sources: 5C4.85 (NGC 4874), 5C4.81 (NGC 4869), 5C4.74, 
5C4.114, 5C4.127, 5C4.42, and 5C4.152. Radio observations 
were performed using the 6 cm and 3.6 cm bands for all of 
the source except 5C4.114. This source being weaker, was 
observed at 20 cm and 6 cm. The sources were observed at 
two frequencies within each band, in order to have 4 frequency 
observations: 4.535 GHz, 4.935 GHz, 8.085 GHz, and 8.465 
GHz. In addition 4.735 GHz observations were performed 
on sources 5C4.85 and 5C4.74. Due to technical issues the 
observing time for the source 5C4.152 was reduced, and it was 
necessary to increase the signal-to-noise ratio in the 3.6 cm 
band. To increase the signal-to-noise ratio, the data at 8.085 
and 8.465 GHz were averaged together and a single image 
at 8.275 GHz was obtained. For the source 5C4.114 we have 
observations at 1.365 GHz, 1.516 GHz, 4.535 GHz and 4.935 
GHz. Details of the observations are reported in Table Q] The 
source 3C286 was used as both primary flux density calibratotQ 
and as absolute reference for the electric vector polarization 
angle. The source 1310+323 was observed as both a phase and 
parallactic angle calibrator. 

We performed standard calibration and imaging using the 
NRAO Astronomical Imaging Processing Systems (AIPS). 
Cycles of phase self-calibration were performed to refine 
antenna phase solutions on target sources, followed by a final 
amplitude and gain self-calibration cycle in order to remove 
minor residual gain variations. Total intensity, I, and Stokes pa- 
rameter Q and U images have been obtained for each frequency 
separately. After cleaning, radio images were restored with a 
naturally weighted beam. The final images were then convolved 
with a Gaussian beam having FWHM = 1.5"xl.5" (~ 0.7x0.7 
kpc). Polarization intensity P - ^U 2 + Q 2 , Polarization angle 
*F = jatan(U, Q) and fractional polarization FPOL - j images 
were obtained from the I, Q and U images. Polarization intensity 
images have been corrected for a positive bias. The calibration 
errors on the measured fluxes are ~ 5%. 



1 we refer to the flux density scale by Baars et al. (1990) 
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Table 1. VLA observations of radio galaxies in the Coma cluster field. 
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Source 


RA 


DEC 


V 


Bandwidth 


Config. 


Date 


Time on Source 




rT?00(Tl 

\J L,\J\J\J ) 




(GHz) 


(MHz) 








5C4.85 


12 59 35.3 


+21 57 36 


8.085 - 8.465 


50 


B 


Jul 06 


2.6 








4.535 - 4.935 


50 


B 


Oct 07 


2.6 








4.735 


50 


B 


May 09 


2.6 


5C4.81 


1 2 59 22 8 


+27 54 40 


8 085 - 8 465 


50 


B 


Jul 06 


2.6 








4.535 - 4.935 


50 


B 


Jul 06 


2.7 








8.465 


50 


C 


Nov 90 


3.8 


5C4.74 


12 58 59.4 


+27 46 46 


8.085 - 8.465 


50 


B 


Jul 06 


2.7 








8.085 - 8.465 


50 


C 


Apr 08 


5.4 








A - A Q^S 


50 


g 


Jul 06 


2.7 








4.735 


50 


B 


May 09 


2.6 


^CA 1 1 A 


i ^ on so & 


_i_98 08 m 
TiO uo vto 


1 .JUJ - 1 . J 1 J 


95 


A 


rw 08 










4.535 - 4.935 


50 


B 


Apr 09 


5.1 


5C4.127 


13 01 20.1 


+28 21 38 


8.085 - 8.465 


50 


B 


Jul 06 


2.6 








4.535 - 4.935 


50 


B 


Oct 07 


2.9 


5C4.42 


12 56 38.2 


+28 16 43 


8.085 - 8.465 


50 


B 


Jul 06 


2.5 








4.535 - 4.935 


50 


B 


Jul 06 


2.6 


5C4.152 


13 03 14.4 


+27 30 06 


8.085 - 8.465 


50 


B 


Aug 06 


1.2 








4.535 - 4.935 


50 


B 


Aug 06 


1.5 



Col. 1: Source name; Col. 2, Col. 3: Pointing position (RA, DEC); Col. 4: Observing frequency; 
Col 5: Observing bandwidth; Col. 6: VLA configuration; Col. 7: Dates of observation; 
Col. 8: Time on source (flags taken into account). 



Table 2. Total and polarization intensity radio images. Images are restored with a beam of 1.5"xl.5" 



Source name 


V 


o-(I) 


o-(Q) 


o-(U) 


Peak brightness 


Flux density 


Pol. flux 




(GHz) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy/beam) 


(mJy) 


(mJy) 


5C4.85 


4.535 


0.018 


0.015 


0.015 


5.8 


75.7 


7.5 




4.735 


0.017 


0.016 


0.016 


5.7 


75.6 


7.6 




4.935 


0.019 


0.017 


0.017 


5.8 


75.2 


7.5 




8.085 


0.032 


0.021 


0.023 


3.3 


38.6 


4.2 




8.465 


0.034 


0.022 


0.022 


3.1 


36.4 


4.0 


5C4.81 


4.535 


0.025 


0.022 


0.024 


2.7 


66.7 


12.1 




4.935 


0.025 


0.022 


0.022 


2.7 


58.2 


10.5 




8.085 


0.022 


0.019 


0.021 


2.4 


25.2 


5.5 




8.465 


0.015 


0.014 


0.014 


2.3 


54.0* 


11.3 


5C4.74 


4.535 


0.025 


0.023 


0.023 


2.7 


17.0 


3.4 




4.735 


0.019 


0.016 


0.016 


2.3 


13.4 


2.7 




4.935 


0.023 


0.021 


0.021 


2.6 


15.6 


3.1 




8.085 


0.015 


0.013 


0.013 


1.7 


9.1 


2.2 




8.465 


0.014 


0.013 


0.013 


1.6 


8.6 


2.1 


5C4.114 


1.365 


0.040 


0.027 


0.027 


11.3 


47.0 


5.9 




1.516 


0.034 


0.021 


0.022 


10.4 


42.6 


5.5 




4.535 


0.014 


0.013 


0.013 


4.2 


16.4 


3.0 




4.935 


0.014 


0.013 


0.013 


3.8 


14.9 


2.8 


5C4.127 


4.535 


0.028 


0.019 


0.017 


56.3 


72.5 


3.9 




4.935 


0.024 


0.020 


0.017 


56.4 


72.9 


3.9 




8.085 


0.023 


0.021 


0.022 


52.2 


62.3 


3.4 




8.465 


0.025 


0.021 


0.023 


51.6 


61.4 


3.3 


5C4.42 


4.535 


0.023 


0.022 


0.022 


7.2 


63.2 


8.0 




4.935 


0.023 


0.022 


0.022 


6.5 


57.8 


7.3 




8.085 


0.022 


0.020 


0.020 


3.8 


33.2 


4.1 




8.465 


0.021 


0.020 


0.019 


3.6 


31.0 


4.0 


5C4.152 


4.535 


0.026 


0.024 


0.025 


7.0 


22.7 


2.9 




4.935 


0.026 


0.025 


0.025 


6.4 


20.3 


2.4 




8.275 


0.021 


0.019 


0.019 


4.1 


12.3 


1.8 



Col. 1: Source name; Col. 2: Observation frequency; Col. 3, 4, 5: RMS noise of the I, Q, U images; 
Col. 7: Peak brightness; Col. 8: Flux density; Col. 9: Polarized flux density. 

"The higher flux measured at 8.465 GHz is derived by combining together B and C array observations. 
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30:00.0 



TO 
£ 

~o 
a> 
Q 



28:00:00.0 



30:00.0 



27:00:00.0 

04:00.0 02:00.0 13:00:00.0 58:00.0 12:56:00.0 

Right ascension 

Fig. 1. Colors: Coma X-ray emission from the ROSAT All Sky Survey in the energy band [0.1, 2.4] kev. Contours: Coma radio 
emission at 1.4 GHz from the NVSS. The beam FWHM is 45"x45", contours start from 1.5 mJy/beam and are spaced by a factor 
of 2. The observed sources are labelled. 



3.2. Radio properties of the observed sources 

In this section the radio properties of the observed sources are 
briefly presented. Further details are given in Table|2] 
Redshift information is available for three out of the seven ob- 
served radio sources. Two of them (5C4.85 and 5C4.8 1) are well 
studied Coma cluster members, while the third one (5C4.127) 
is associated with a background source. Although the redshift 
is not known for the other four radio sources, they have not 
been identified with any cluster member down to very faint op- 
tical magnitudes: M r > -15 (see Miller et al. 2009). This indi- 
cates that they are background radio sources, seen in projection 
through the cluster. In the following the radio emission arising 
from the selected sample of sources is described together with 
their main polarization properties. In the fractional polarization 
images (from Fig.|2]to[8]) pixels with errors larger than 10% were 
blanked. 

5C4.85 - NGC 4874 

This a cluster source, optically identified with the Coma cen- 
tral cD galaxy NGC 4874 (see e.g. Mehlert et al. 2000). It is a 
Wide Angle Tail radio galaxy, whose maximum angular exten- 
sion is ~ 30", corresponding to ~ 15 kpc. The angular extension 
of the two lobes individually is larger at the lowest frequency. 
The northern lobe shows a mean fractional polarization of 10% 



and 1 1 % at 4.535 and 8.465 GHz respectively, while the western 
lobe is less polarized ( ~7% at both frequencies). In Fig. [2] the 
radio emission is shown at 4.535 and 8.465 GHz. 

5C4.81 - NGC 4869 

This source has been studied in detail by Dallacasa et al. (1989) 
and Feretti et al. (1995). It is associated with the giant elliptical 
galaxy NGC4869. 5C4.81 has a Narrow Angle Tail radio mor- 
phology, and its angular size in the images obtained here is 55" 
(25 kpc). The mean fractional polarization in the tail is 18% at 
4.535 GHz and 21% at 8.465 GHz. In Fig.[3]the radio emission 
is shown at 4.535 and 8.465 GHz. 

5C4.74 

The source 5C4.74 consists of 5C4.74a and 5C4.74b, the two ra- 
dio lobes of a FRII radio source. Its redshift is unknown, and no 
optical identification has been found, either with a Coma cluster 
member (Miller et al. 2009) nor with a background radio source. 
From this we conclude that it is a distant background source. 
The northeastern lobe has a fractional polarization of ~28% and 
35% respectively at 4.535 and 8.465 GHz, while the southwest- 
ern lobe is less polarized (~ 19% at 4.535 GHz and ~ 23% at 
8.465 GHz). 

In Fig.|4]the radio emission is shown at 4.535 and 8.465 GHz. 
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Fig. 2. Source 5C4.85. Total intensity radio contours and polarization vectors at 4.535 GHz (left) and 8.465 GHz (right). The bottom 
contour corresponds to a 3cr noise level, contours are then spaced by a factor of 2. E vectors are superimposed: the orientation 
indicates the direction of the E field, while the line length is proportional to the fractional polarization intensity (1" corresponding 
to 10%). 




K.5 2i0 Zl{ Hfl 21.5 21j0 2Cl5 3ClO 1M 2i0 Zlf 22j0 21.5 21.0 30l5 3Cl0 

HGHT ASCENSION (J200O| HG HT ASC EN30 N (J200O) 



Fig. 3. Source 5C4.8 1 . Total intensity radio contours and polarization vectors at 4.535 GHz (left) and 8.465 GHz (right). The bottom 
contour corresponds to a 3cr noise level, and contours are then spaced by a factor of 2. E vectors are superimposed: orientation 
indicates the direction of the E field, while line length is proportional to the fractional polarization intensity (1" corresponding to 
10%). 



5C4.114 

5C4.114 is a FRI radio source, with angular size of ~ 15". Its 
redshift is unknown, and no optical identification either with a 
Coma cluster galaxy (Miller et al. 2009) nor with a background 
galaxy has been found, indicating that 5C4.114 has a redshift 
greater than 0.023. The southern lobe appears brighter than the 
northern one. The source fractional polarization is ~ 13% at 
1.365 GHz and ~ 19% at 4.935 GHz. 

In Fig.[5]the radio emission is shown at 1.365 and 4.935 GHz. 
5C4.127 

5C4.127 is a QSO located at z= 1.374 (Veron-Cetty & Veron, 
2001). Observations presented here show that in addition to a 
bright nucleus the source has a weak extension in the E-W di- 
rection of ~ 16 " (~ 136 kpc) at both of the observing frequency 
bands. The extended component has a mean fractional polariza- 
tion of 13% at 4.535 GHz and 14% at 8.465 GHz, while the 



nucleus is polarized at the 3% level. In Fig.|6]radio contours of 
the source and polarization vector images are shown. 

5C4.42 

5C4.42 is a FRII-type radio source. Redshift information is 
not available in the literature and no optical identification has 
been found. The same arguments explained above for the source 
5C4.74 let us conclude that it is a background radio source. The 
source is composed by a weakly polarized core and two lobes 
that extend for ~ 25" in the southwest and northeast directions. 
The lobes show a mean fractional polarization of ~ 13% at both 
4.535 GHz and 8.465 GHz. In Fig.Qradio contours and vector 
polarization images of the source are shown. 

5C4.152 

5C4.152 is a FRII type Radio Galaxy. No redshift is available 
in the literature for this source. The same arguments explained 
above for the source 5C4.74 let us conclude that it is a back- 
ground radio source. It is composed of a core having a fractional 
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125902 01 CO 5359 53 57 12 59 OS 01 00 5359 53 57 

HGHTAECENaON(JBMO) HG KT AECE NEIO N &I3D00) 



Fig. 4. Source 5C4.74. Total intensity radio contours and polarization vectors at 4.535 GHz (left) and 8.465 GHz (right). The bottom 
contour corresponds to a 3cr noise level, contours are then spaced by a factor of 2. E vectors are superimposed: the orientation 
indicates the direction of the E field, while the line length is proportional to the fractional polarization intensity (1" corresponding 
to 10%). 




Fig. 5. Source 5C4.114. Total intensity radio contours and polarization vectors at 1.365 GHz (left) and 4.935 GHz (right). The 
bottom contour corresponds to a 3cr noise level, contours are then spaced by a factor of 2. E vectors are superimposed: the orientation 
indicates the direction of the E field, while the line length is proportional to the fractional polarization intensity (1" corresponding 
to 10%). 



polarization of a few percent and two lobes that extend for ~ 
28" north-south. The lobes show a mean fractional polarization 
of ~ ~ 13% 4.535 GHz and 15% at 8.275 GHz. In Fig.[8]radio 
contours and vector polarization images of the source are shown. 

4. RM: fits and errors 

In this section the procedure used to derive the RM from 
the radio observations is explained. We used the PACERMAN 
algorithm (Polarization Angle CorrEcting Rotation Measure 
ANalysis) developed by Dolag et al. (2005). The algorithm 
solves the nn ambiguity in low signal-to-noise regions by ex- 
ploiting the information of nearby reference pixels, under the 
assumption that the reference pixel is connected to the nearby 
areas as far as the polarization angle gradient is below a certain 
threshold in all of the observed frequency maps simultaneously. 



We considered as reference pixel those with a polarization angle 
uncertainty less than 7 degrees, and fixed the gradient threshold 
to 15 degrees. An error of 7 degrees in the polarization angle 
corresponds to 3<x level in both U and Q polarization maps si- 
multaneously. We allowed PACERMAN to perform the RM fit 
if at least in 3 frequency maps the above mentioned conditions 
were satisfied. The resulting RM images are shown in Fig. [9] 
[10l[ni[T2l[T3][T4]and[T5]overlaid on the total intensity contours 
at 4.935 GHz. In the same figures we also provide the RM dis- 
tribution histograms and the RM fits for selected pixels in the 
map. The linear trend of *F versus A 2 and the good fits obtained 
clearly indicate that the Faraday rotation is occurring in an ex- 
ternal screen. 

From the RM images we computed the RM mean {{RM)) and 
its dispersion {(jRM,obs)- There are two different types of errors 
that we have to account for: the statistical error and the fit er- 
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Fig. 7. Source 5C4.42. Total intensity radio contours and polarization vectors at 4.535 GHz (left) and 8.465 GHz (right). The bottom 
contour corresponds to a 3cr noise level, contours are then spaced by a factor of 2. E vectors are superimposed: the orientation 
indicates the direction of the E field, while the line length is proportional to the fractional polarization intensity (1" corresponding 
to 10%). 



ror. The statistical errors for (RM) and for cr RMobs is given by 
o~RM,obsl y/ftb and cr RMo i,sl y/2rib respectively, where Hb is the 
number of beams over which the RM has been computed. The 
statistical error is the dominant one, while the error of the fit has 
the effect of increasing the real (Trm- Thus, in order to recover 
the real standard deviation of t he observed RM distributio n we 
have computed the cr RM dec as -yJcr RMobs — Median(Errji,) 2 . with 

Median(Err f lt ) being the median of the error distribution. The 
fit error has been estimated with Monte Carlo simulations. We 
have extracted n R values, from a random Gaussian distribution 
having <r = cr RMo i, s and mean = (RM), we have then added to 
the extracted values a Gaussian noise having cr„ i se = Err/i t , 
in order to mimic the effect of the noise in the observed RM 
images. We have computed the mean and the dispersion (cr J(m ) 
of these simulated quantities and then subtracted the noise from 

the dispersion obtaining cr sim j ec = ^j(T 2 jm - cr 2 nohe . We have thus 

obtained a distribution of a- sim j ec and means. The standard de- 
viation of the (Tsim,dec distribution is then the fit error on cr RMc i ec 



while the standard deviation of the mean distribution is the fit 
error on (RM). We checked that the mean of both distributions 
recover the corresponding observed values. In Table[3]we report 
the RM mean, the observed RM dispersion (o- RMo i, s ), the value 
of cr R M,dec (hereafter simply cr RM ), with the respective errors, the 
average fit error (Errf,,), and the number of beam over which the 
RM statistic is computed («/,). 



4.0.1. The source 5C4.74 

The value of (RM) that we have derived for the source 5C4.74 is 
quite high compared with the values found for the other sources 
in this cluster and it is also higher than the values obtained in 
other clusters for sources at similar distances from the cluster 
center (e.g. Clarke et al. 2004). The level of polarization of this 
source is also quite high compared to the other sources (see Sec. 
I3.2I ). We note its position southwest of the cluster core, in the 
direction of the sub-group NGC4839 that is currently merging 
with the Coma cluster (Feretti & Neumann 2006). One possi- 
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Fig. 8. Source 5C4.152. Total intensity radio contours and polarization vectors at 4.535 GHz (left) and 8.465 GHz (right). The 
bottom contour corresponds to a 3<x noise level, contours are then spaced by a factor of 2. E vectors are superimposed: the orientation 
indicates the direction of the E field, while the line length is proportional to the fractional polarization intensity (1" corresponding 
to 10%). 



bility is that the magnetic field and/or the thermal gas has been 
compressed and ordered in this region, increasing the observed 
polarization flux and making (RM) peculiar in this position. This 
might imply that more sophisticated models, that include devi- 
ations from a spherical symmetry could give a better represen- 
tation of the gas density profile. We note however that the X- 
ray analysis performed in the literature by Briel et al. (1992) 
shows that the spherical /?-model is a good representation of the 
cluster X-ray surface brightness, indicating that deviations from 
spherical symmetry are small. The cr RM value is fully compatible 
with the trend suggested by the other sources. This could be ex- 
plained in the proposed scenario, if the magnetic field has been 
compressed and ordered. In the following analysis we will use 
mainly the cr RM to infer the magnetic field strength and structure. 
Nonetheless, results will also be presented excluding this source 
from our analysis. 

4.1. Galactic contribution 

The contribution to the Faraday RM from our Galaxy may 
introduce an offset in the Faraday rotation that must be removed. 
This contribution depends on the galactic positions of the 
observed sources. The Coma cluster Galactic coordinates are 
/ = 58° and b = 88°. The cluster is close to the galactic north 
pole, so that Galactic contribution to the observed RM is likely 
negligible. However, in order to estimate this contribution the 
average RM for extragalactic sources located in projection 
nearby the Coma cluster region has been computed using the 
catalogue by Simard-Normadin et al. (1981). The sources have 
been weighted by the inverse of the distance from the Coma 
cluster center, and the results in a region of 25x25 degrees 2 
centered on the cluster gives a Galactic contribution of -0.15 
rad/m 2 . This small contribution is thus completely negligible 
and has been ignored in the following analysis. 



4.2. RM local contribution 

We discuss here the possibility that the RM observed in radio 
galaxies are not associated with the foreground ICM but may 



arise locally to the radio source (Bicknell et al. 1990, Rudnick & 
Blundell 2003), either in a thin layer of dense warm gas mixed 
along the edge of the radio emitting plasma, or in its immediate 
surroundings. There are several arguments against this interpre- 
tation: 

- the trend of RM versus the cluster impact parameter in 
both statistical studies and individual cluster investigations 
(Clarke et al. 2001, 2004; Feretti et al. 1999; Govoni et al. 
2005); 

- the Laing-Garrington effect (Laing 1988, Garrington et al. 
1988, Garrington & Conway 1991). This effect consists of 
an asymmetry in the polarization properties of the lobes 
of bright radio sources with one-sided, large scale jets. 
The lobe associated with the jet that is beamed toward the 
observer is more polarized than the one associated with the 
counter-jet that points away from the observer. This effect 
can be explained if we assume that the radio emission from 
the two lobes cross different distances through the ICM, and 
therefore the emission from the counter-lobe is seen through 
a greater Faraday depth, causing greater depolarization. This 
means also that the observed polarization properties of the 
source are strongly influenced by the ICM. 

- statistical tests on the scatter plot of RM versus polarization 
angle for the radio galaxy PKS 1246-410 (Ensslin et al. 
2003); 

- the relation between the RM and the cooling flow rate in 
relaxed clusters (Taylor et al. 2002). 

This observational evidence allows us to conclude that the main 
contribution to the RM occurs in a Faraday screen located out- 
side the radio sources. The only contribution that could contam- 
inate the cluster Faraday screen is the contribution from the dis- 
turbed ISM in the central parsecs of the host elliptical galaxy 
where the inner jet has been found to have RMs up to thousands 
of radians per square meter (Zavala & Taylor 2004). We removed 
the core from the RM analysis in order to avoid any contribution 
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Fig. 9. 5C4.85: Top left: The RM fit is shown in color along with total intensity radio contours 4.935 GHz. The bottom contour 
correspond to the 3<x noise level and contours are then spaced by a factor of l.Top right distribution histogram of the RM values. 
Bottom: fits of polarization angle versus A 1 in four representative pixels. 



Table 3. Rotation Measures values of the observed sources 



Source 


Projected distance 


n. of beams 


(RM) 


O" RM,obs 


Err /ft 


CRM 




kpc 




rad/m 2 


rad/m 2 


rad/m 2 


rad/m 2 


5C4.85 


51 


35 


-256±50 


303 


46 


299±36 


5C4.81 


124 


56 


-120±22 


166 


48 


159+17 


5C4.74 


372 


10 


372±51 


154 


44 


148±41 


5C4.114 


532 


16 


51±4 


16 


2 


16±3 


5C4.127 


919 


7 


21±30 


65 


36 


54±26 


5C4.42 


1250 


33 


6±12 


56 


43 


36+11 


5C4.152 


1489 


4 


32±27 


37 


28 


24±21 



Col. 1 : Source name Col. 2: Source projected distance from the X-ray cluster center; 
Col. 3: number of beams over which RMs are computed; 
Col. 4: Mean value of the observed RM distribution; 
Col. 5: Dispersion of the observed RM distribution; 

Col. 6: Median of the RM fit error; Col 7: Dispersion of the RM distribution after noise deconvolution. 



of this kind (out to a distance of ~5" from the core). 
The ICM origin of the observed RM is also confirmed by the data 
presented here (Tab.0: the trend of <x SM exhibits a decrease with 
increasing cluster impact parameter. Values of (RM) £0 and dif- 
ferent among sources located at different projected distances to 
the cluster center indicate that the magnetic field substantially 
changes on scales larger than the source size, while small RM 
fluctuation can be explained by magnetic field fluctuation on 
scales smaller than the source size. Thus in order to interpret 
correctly the RM data we have to take into account magnetic 
field fluctuations over a range of spatial scales, i.e., we have to 
model the magnetic field power spectrum. 



5. The magnetic field model 

In this section the magnetic field model adopted to simulate 
RM images is described. We used the approach suggested by 
Murgia et al. (2004), where the magnetic field is modelled as 
a 3-dimensional multi-scale model and its intensity scales with 
radius, following the gas distribution. 

5.1. The magnetic field power spectrum 

In order to study the magnetic field of the Coma cluster, we con- 
sidered a 3D vectorial magnetic field model. 
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Fig. 10. 5C4.81:7bp left: The RM image is shown in color along with total intensity radio contours at 4.935 GHz. Contours start at 
3<x and increase by factors of 2. Top right distribution histogram of the RM values. Bottom: fits of polarization angle versus A 2 in 
four representative pixels. 



Simulations start considering a power spectrum for the vector 
potential A in the Fourier domain: 

\A k \ 2 K (4) 

and extract random values of its amplitude A and phase (p. A 
is randomly extracted from a Rayleigh distribution (in order to 
obtain a Gaussian distribution for the real magnetic field com- 
ponents, as frequently observed), while <p varies randomly from 
to 2jt. The magnetic field components in the Fourier space are 
then obtained by: 

B(k) = ik x A(k). (5) 

The field components B, in the real space are then derived us- 
ing 3D Fast Fourier Transform. The resulting magnetic field is a 
multi-scale model with the following properties: 
1) 

V • B = 0, 

2) the magnetic field energy density associated with each com- 
ponent Bit is: 

\B k \ 2 = C\k-\ 

n = f — 2, where C 2 is the power spectrum normalization, 

3) Bj has a Gaussian distribution, with (B,) = 0, cr B = (B 2 ), 

4) B has a Maxwellian distribution. 



The magnetic field model scales with the radius as described in 
Sec.EH 

We define A = 25 as the physical scale of the magnetic field fluc- 
tuations in the real space. Thus in order to determine the mag- 
netic field power spectrum in the cluster, we have to determine 
three parameters: A„„„, A max and n. It is worth noting that a de- 
generacy arises between A max and n (the higher n is, the lower 
A max is, in order to produce the same RM, see also Sec. 17. U . 

5.2. The magnetic field radial profile 

There are several indications that the magnetic field intensity de- 
creases going from the center to the periphery of a cluster. This 
is expected by magneto-hydrodynamical simulations (see e.g. 
Dolag et al. 2008) and by spatial correlations found in some clus- 
ters between thermal and non-thermal energy densities (Govoni 
etal.2001). 

We assume that the cluster magnetic field follows the thermal 
component radial distribution according to: 

<B>(r) = <B„>(^)" (6) 

where (Bo) is the mean magnetic field strength at the cluster cen- 
ter. 

In order to obtain the desired magnetic field radial profile we 
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Fig. 11. 5C4.74:Top left: The RM image is shown in color along with total intensity radio contours at 4.935 GHz. Contours start at 
3<x and increase by factors of 2. Top right distribution histogram of the RM values. Bottom: fits of polarization angle versus A 2 in 
four representative pixels. 



have operated directly in the real space. Strictly, this operation 
should be performed in the Fourier space, by convolving the 
spectral potential components with the shaping profile. It has 
been proved that these two approaches give negligible differ- 
ences (see Murgia et al. 2004). 

When the magnetic field profile is considered, two more param- 
eters have to be determined: 77 and (Bo). These two parameters 
are degenerate (the higher is (Bo) the higher is rj). In fact higher 
central value of the field require a steeper spectrum in order to 
produce the same RM. 

The adopted magnetic field model has then a total of 5 free 
parameters: A„„„, K max , n, rj and (-Bo), and is subject to two de- 
generacies: A max - n and rj and (Bq). 

Fitting all of these five parameters simultaneously would be the 
best way to proceed, but it is not feasible here, due to the com- 
putational burden caused by the Fourier Transform inversion. 
Indeed we have to simulate a large volume ~ 3 3 Mpc 3 with a 
sub-kiloparsec pixel-size. 

The aim of this work is to constrain the magnetic field radial 
profile, and for this reason the sources were selected in order 
to sample different regions at different impact parameters. This 
allows us to reach a good sensitivity to the RM at different dis- 
tances from the cluster center. 

We proceed as follows: we perform 2D simulations with differ- 
ent magnetic field power spectra in order to recover the RM sta- 



tistical indicator that are sensitive to the magnetic field power 
spectrum (Sec. 17. U . From this analysis we derive the power 
spectrum that best reproduces the observations. We then perform 
3D magnetic field simulations varying the values of Bq and 77 and 
derive the magnetic field profile that best reproduces the RM ob- 
servations (Sec. 17.3b . 



6. Comparing observations and simulations 

A tricky point when observations and simulations are compared 
is the correct evaluations of the errors and uncertainties that this 
process is subject to. The simulations we present in this paper 
start from a random seed and generate 2D and 3D magnetic 
fields. From these fields simulated RM images are obtained, and 
then compared with those observed in order to constrain the 
magnetic field properties. It is worth noting that for a given mag- 
netic field model, the RM in a given position of the cluster varies 
depending on the initial seed of the simulation, so that different 
realizations of the same model will correspond to different val- 
ues of (RM) and cr RM at that position. Given the limited area 
covered by RM observations, the random nature of the field can- 
not be neglected in our analysis. 

We adopt the following approach to compare observations and 
simulations: once the simulated RM image is obtained for a 
source, it is convolved with a Gaussian function having FWHM 
equal to the beam FWHM of the observed image. The simu- 
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Fig. 12. 5C4.114:7bp /e/f: The RM image is shown in colors along with total intensity radio contours at 4.935 GHz. Contours start 
at 3<x and increase by factors of 2. Top right distribution histogram of the RM values. Bottom: fits of polarization angle versus A 2 in 
four representative pixels. 



lated RM image is then blanked in the same way as the observed 
RM image. This ensures that simulations are subject to the same 
sampling bias that we have to deal with when obtaining the RM 
from observations. The comparison between the observed RM 
images and those simulated is performed with the x 1 distribu- 
tion, by computing: 



X 



-S 



1=1 



Errj, 



(7) 



where i indicates the source, C t s refers to a generic observed 
quantity, while <C, iS ,' m ) is the same quantity averaged over the 
different equivalent numerical simulation that start with different 
random seeds, and Errl, refers to the error of C b s - 



7. Determining the magnetic field from RM 
observations 

Here we describe how the magnetic field power spectrum has 
been investigated. 

7.1. Constraining the magnetic field power spectrum 

Several observational quantities can be useful to constrain some 
properties of the magnetic field power spectrum. In particular: 



- Both (RM) and cr RM scale linearly with the magnetic field 
strength, while they have different trends with n and A max , 
which are degenerate parameters. The ratio \(RM)\/ctrm can 
thus be used to investigate the magnetic field power spectrum 
(see also Fig. 3 in Murgia et al. 2004). 

- The minimum scale of the magnetic field fluctuation, A m ,„, 
affects the depolarization ratio (DP ratio) at two different 



frequencies (i.e. DP - 



) and the o-rm- Both DP and 



o-RM are in fact determined by the magnetic power on the 
small spatial scales. This parameter can be thus be derived 
by studying high resolution polarization images. 
- It has been demonstrated that the magnetic field auto- 
correlation function is proportional to the RM auto- 
correlation function (Ensslin & Vogt 2003). Since the power 
spectrum is the Fourier transform of the auto-correlation 
function, it is possible to study the 3D magnetic field power 
spectrum starting from the power spectrum of the RM im- 
ages. 

We simulated 2D magnetic field models with different power 
spectra and compared simulated RM images and DP with the 
corresponding observable quantities. In these simulations the 
power spectrum normalization is set independently for each 
sourc^l The computational grid is 512 x 512 pixel 2 and the 



2 In these 2D RM simulations the radial profile of the magnetic field 
is not accounted for. This implicitly assumes that the mean magnetic 
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Fig. 13. 5C4.127 :Top left: The RM image is shown in color along with total intensity radio contours at 4.935 GHz. Contours start 
at 3cr and increase by factors of 2. The circle indicates the area masked in the RM analysis. Top right distribution histogram of the 
RM values. Bottom: fits of polarization angle versus A 2 in four representative pixels. 



pixel-size was fixed to 0.2 kpc. This guarantees that each beam 
is represented by three pixels in the grid. The resulting field of 
view is then ~ 100x100 kpc 2 , that is enough to recover the pro- 
jected size of the sources and to properly sample the large power 
spectrum scales. Simulations were convolved with a Gaussian 
function having FMHM equal to the beam of the observations. 

7.1 .1 . The A max -n plane 

In order to illustrate the degeneracy existing between A max and 
n, 2D magnetic field models with different power spectra have 
been simulated. We allowed the parameter n to vary in the range 
[0.5; 4.5] and A max in the range [5; 600] kpc, with steps of 0.06 
and 9.3 kpc respectively. We derived simulated RM images for 
each combination of these parameters and calculated the RM 
ratio as: 

\{RM)\ 

RMratio = - (8) 

CRM 

in a region of 15 x 15 kpc 2 , comparable to the regions where 
RM has been observed. In Fig.[l6]the values of the RM ratio for 
the simulated RM images is shown in colors, as a function of 
both n and A max . The RM ratio was computed for the observed 



field strength is not dramatically varying over the scale of the source. 
This is a reasonable assumption since the linear sizes of the sources at 
the Coma redshift are in fact much smaller that the cluster core radius. 



source that have a RM signal-to noise ratio >3 both for cr RM and 
(RM), i.e for the sources 5C4.85, 5C4.81, 5C4.74 and 5C4.1 14. 
The resulting values are shown in Fig. [16] (black line). 
The plot shown in Fig [16] shows what A max - n degeneracy 
means: the same value of the RM ratio can be explained with 
different power spectra. There are, as expected, two asymptotic 
trends. In fact, if the magnetic field power spectrum is flat (e. 
g. n <2), the bulk of the magnetic field energy is on the small 
scales, and thus the effect of increasing A max is negligible after 
a certain threshold, that in this case is achieved for A max ~300 
kpc. As the power spectrum steepens the bulk of the energy 
moves to large scales, and thus as A max increases, the energy 
content also increases sharply. This is the reason of the second 
asymptotic trend that is shown in the plot: as n increases A max 
decreases faster and faster. As n approaches the value of ~ 1 1/3 
(Kolmogorov power spectrum), the observed data constrain 
A max to be ~20 - 40 kpc. 



7.2. Structure function, auto-correlation function and 
Multi-Scale-Statistic 

In order to constrain more precisely the estimate of the magnetic 
field power spectrum parameters indicated by the previous anal- 
ysis we have investigated the statistical properties of the RM im- 
ages individually. We have fixed n = 11/3, corresponding to the 
Kolmogorov power law for turbulent fields. This choice is mo- 
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Fig. 14. 5C4.42:Top left: The RM image is shown in color along with total intensity radio contours at 4.935 GHz. Contours start at 
3<x and increase by factors of 2. Top right distribution histogram of the RM values. Bottom: fits of polarization angle versus A 2 in 
four representative pixels. 



tivated by both observational and theoretical works. Schuecker 
et al. (2004) analyzed spatially-resolved gas pseudo-pressure 
maps of the Coma galaxy cluster deriving that pressure fluctua- 
tions in the cluster center are consistent with a Kolmogorov-like 
power spectrum. Furthermore, cosmological numerical simula- 
tions have recently demonstrated that 3D power spectrum of the 
velocity field is well described by a single power law out to at 
least one virial radius, with a slope very close to the Kolmogorov 
power law (Vazza et al. 2009a, b). 

The range of values of A max is suggested by the previous anal- 
ysis (see Fig.fTol). In order to choose the best parameters in that 
range, and to find the best value for A,„,„, we simulated RM im- 
ages and used two different statistical methods to compare the 
observed RM images to the simulated ones: 

1 . We calculated the auto-correlation function and the structure 
function of the observed RM image, and then compared them 
with the simulated RM image. The RM structure function is 
defined as follows. 

S (dx, dy) = ([RM(x, y) - RM(x + dx,y + dy)]\ xj) , (9) 

where = Ofey) indicates that the average is taken over all the 
positions (x,y) in the RM image. Blank pixels were not con- 
sidered in the statistics. The structure function S(r) is then 
computed by radially averaging S (dx, dy) over regions of in- 
creasing size of radius r = yjdx 1 + dy 2 . S(r) is thus sen- 



sitive to the observable quantity cr RM over different scales. 
The auto-correlation function is defined as: 

A(dx, dy) = ( [RM(x, y)RM(x + dx,y + dy)]) {xj) ( 1 0) 

Since A(0) = {RM 2 } - cr 2 RM + {RM) 2 , the auto-correlation 
function is sensitive to both {RM) and the ctrm- 
2. We computed a Multi-Scale Statistic, namely we computed 
{RM) and cr RM over regions of increasing size in the ob- 
served RM images and compared them with the same values 
obtained in the simulated images. The smallest region over 
which {RM) and cr RM are computed corresponds to a box of 
0.4 x 0.4 kpc size. The box side is then increased by a factor 
two until the full source size is reached. We note that this 
approach is sensitive to both {RM) and cr RM over different 
spatial scales, and is thus a useful tool to discriminate among 
different power spectra. This indicator differs from the S (r) 
and A(r) in that as r increases, the number of pixels useful 
for computing the Multi-Scale Statistic increases, giving a 
robust statistical estimate on large scales. 

For each source we simulated different power spectra varying 
A m ;„ from 1 kpc to 5 kpc and A max from 20 to 40 kpc. For every 
power spectrum and for each source we realized ten different 2D 
RM images, and compared the statistics (S(r), A(r) and Multi- 
Scale Statistic) with the observed ones by computing Eq.|7] This 
approach makes it possible to discriminate the best power spec- 
trum model compatible with our data. In this 2D analysis we 
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Fig. 15. 5C4.152:rop left: The RM image is shown in color along with total intensity radio contours at 4.935 GHz. Contours start 
at 3<x and increase by factors of 2. Top right distribution histogram of the RM values. Bottom: fits of polarization angle versus A 2 in 
four representative pixels. 



focused on the sources 5C4.85, 5C4.81, 5C4.74 and 5C4.114, 
whose RM images have signal-to-noise ratio >3. Each source 
was fitted separately. The total x 2 was then computed by sum- 
ming the individual values obtained for each source. The mini- 
mum value of the total x Z ,,t corresponds to the power spectrum 
model characterized by A„ MA = 34 kpc and A,„,„ = 2 kpc. We 
show in Figs. [17] and [18] the Kolmogorov power spectrum model 
that best reproduces the observed RM images. In Appendix lAl 
similar Figures obtained with other power spectrum models are 
shown. 

7.2.1 . A„„„ and Fractional polarization 

It has been demonstrated (Burn 1966, see also Laing 2008) that 
when Faraday Rotation occurs the fractional polarization FPOL 
can be related to the fourth power of the observing wavelength A 
according to the Burn law: 

p 

FPOL = -i oc exp(-kA A ). (11) 

SinceF POL is sensitive to the minimum scale of the power spec- 
trum, A„„„, Eq.QT]can be used to constrain it. We fitted Eq.fTTIto 
our observations and to our simulations, performed with differ- 
ent values of A m ,„. These fits indicate that a Kolmogorov power 
spectrum with the best agreement with observations is achieved 
for A mm ~ 2 kpc, confirming the result from the previous anal- 
ysis. As an example, we show in Fig.[19]these fits performed on 



the source 5C4.114, where the effect of depolarization is more 
appreciable, thanks to the 20 cm observations. 

7.3. The magnetic field profile 

The results obtained from the previous section indicate the 
power spectrum that is able to best reproduce the observed RM 
images. In order to investigate the magnetic field radial profile 
we simulated 3D Kolmogorov power spectra, with K„ wx - 34 
kpc and A m ,„ = 2 kpc, as derived from the 2D analysis (Sec. 
17. U . A computational grid of 2048 3 pixels was used, and the 
pixel-size was fixed to 0.5 kpc. This guarantees that the Nyquist 
criteria is satisfied for A mln and that fluctuations on scales 
~ A max are also well represented in the cube. For each of these 
simulations Eq. [2] was integrated numerically, with a step of 
0.5 kpc along the line of sight. The limits of the integral in Eq. 
[2] were [0; 10r c ] for the cluster's sources 5C4.85 and 5C4.81 
and [-10r e ; 10r c ] for the other sources in the background 
of the cluster. The simulated field of view covers an area of 
1024x1024 kpc 2 , thus the cube has been replicated to achieve 
a field of view that is large enough to reach the farthest source 
(5C4.152). 

The integration was repeated by varying the parameter Bo in 
the range [0.1; 11] pG, with a step of ~ 0.17 pG, and r\ in the 
range[-0.2; 2.5] with a step of 0.04. For each combination of Bo 
and r] a RM simulated image was thus obtained covering the full 
cluster area. 

We extracted from this RM image seven fields, each lying 
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Fig. 17. Fits to the RM images for the Kolmogorov power spectrum that best reproduces the observed RM (n=l 1/3, A m ,„ =2 kpc, 
A mflA =34 kpc) globally for the sources used in the 2D analysis (see Sec. 17. U . From top to bottom: fit to the <x RM and (RM). 
Red diamonds represent the observed statistic, the cyan line represents the mean taken over ten different realizations of the same 
magnetic field power spectrum, and the blue lines represent the scatter in the simulations. 



in the plane of the sky in the same position of the observed 
sources, and having the same size of the observed RM images. 
The simulated RM images were convolved with a Gaussian 
beam having FWHM=0.7x0.7 kpc, in order to have the same 
resolution of the observations. Finally the simulated RM fields 
were blanked in the same way as the corresponding RM source. 
The result of this integration is, for each combination of (B^rj), 



a set of seven simulated RM images, that are subject to the same 
statistical biases of the observed images. 

This process was repeated 50 times, each starting from a 
different random seed to generate the magnetic field power 
spectrum model. 

For each source and for each value of (Bq;ti) a simulated RM 
image was obtained for every realization of the same power 
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Fig. 18. Fits to the RM images for the Kolmogorov power spectrum that best reproduces the observed RM (n=l 1/3, A m ,„ =2 kpc, 
t^max =34 kpc) globally for the sources used in the 2D analysis (see Sec. 17. 11 1. Fit to the S(r) (top) and A(r) (bottom). Red diamonds 
represent the observed statistic, the cyan line represents the mean taken over ten different realizations of the same magnetic field 
power spectrum, and the blue lines represent the scatter in the simulations. 



spectrum model. The mean and the standard deviation of the 
o~RM,sim(Bo, t]) was computed from the simulated RM images, 
and then the x 2 was obtained (Eq. |7J. The resulting x 2 plane is 
shown in Fig. [20] The minimum value is achieved for Bo =4.7 
fiG and rj = 0.5, but the 1-cr confidence level of the^ 2 indicates 
that values going from Bo =3.9 fiG and r\ = 0.4, to Bo =5.4 
fiG and rj = 0.7, are equally representative of the magnetic 



field profile, according to the degeneracy between the two 
parameters. Magnetic field models with a profile flatter than 
77 < 0.2 and steeper than 77 > 1.0 are excluded at 99% confidence 
level, for any value of (Bo). Also magnetic field models with 
(B ) < 3.1/jG and <B > > 6.5 /jG are excluded at the 99% 
confidence level for any value of 77. It is interesting to note that 
the best models include rj = 0.5, the value expected in the case 
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Fig. 19. Fits to the Burn law. Points refer to observed data, while 
the red line is the fit obtained from observations. Dashed lines 
refer to the fits obtained from three different models, with differ- 
ent values of A ml „, as reported in the bottom left corner of the 
plot. 



of a magnetic field energy density decreasing in proportion 
to the gas energy density (assuming a constant average gas 
temperature), and rj = 0.67, expected in the case of a magnetic 
field frozen into the gas. In the latter case the corresponding 
value of Bq is ~5.2 pG. 

The knowledge of the magnetic field strength and structure 
in the ICM has strong implications for models explaining the 
formation of diffuse radio sources like radio halos. Testing 
the different models proposed in the literature is beyond the 
scope of this paper. We point out, however, that cosmological 
simulations recently performed by Donnert et al. (2009b) have 
shown that it is possible to test a class of these models once the 
magnetic field profile is known. 



7.4. Results excluding the source 5C4.74 

The same procedure described above has been repeated exclud- 
ing the source 5C4.74 (see Sec. 14.0. 11 1. The minimum value for 
the x 2 is again achieved with a model characterized by Bq =4.7 
pG and rj = 0.5. This is not surprising since the fit is computed 
based on <trm, and the source is sampled with only 10 beams. 
In order to investigate possible effects arising from the interac- 
tion of the Coma cluster with the NGC4839 group RM images 
of more sources in this region would be required. 



8. Comparison with other estimates 

In the literature there is a long-standing debate on the magnetic 
field strength derived from the RM analysis compared to the 
equipartition estimate and to the Inverse Compton hard X-ray 



emission. The discrepancy may arise from the different (but not 
incompatible) assumptions, and, moreover, are sensitive to the 
magnetic field on different spatial scales. Assuming the mag- 
netic field models derived in the previous section, it is possible to 
derive an estimate that is comparable with equipartition values, 
and with the Inverse-Compton detection as well as with the up- 
per limits derived from new hard X-ray observations. In order to 
obtain a value that is directly comparable with the equipartition 
magnetic field estimate, we have to derive the average magnetic 
field strength resulting from our RM analysis over the same vol- 
ume assumed in the equipartition analysis, that is ~1 Mpc 3 . The 
magnetic field model resulting from our RM analysis gives an 
average magnetic field strength of ~ 2 pG, consistent with the 
equipartition estimate derived from the radio halo emission ( 0.7 
- 1 .9 pG Thierbach et al. 2003), despite the different assumptions 
that these two methods require. 

The Inverse Compton hard X-ray emission has been observed 
with the Beppo Sax satellite. Its field of view is ~ 1.3°, corre- 
sponding to ~ 2.2 x 2.2 Mpc 2 at the Coma redshift. We com- 
puted the average value of the magnetic field over the same vol- 
ume sampled by Beppo Sax. We obtained ~0.75 pG when the 
best model is assumed, that is a factor four higher than the value 
derived from Hard-X ray observations (Fusco Femiano et al. 
2004). We note however that models compatible with our data 
within 1-cr of the x 1 give values slightly different, going from 
0.9 to 0.5 pG. The steepest magnetic field model that is com- 
patible with our data at 99% confidence level (Bo ~ 6.4 pG, 
rj = 0.95) gives 0.2 pG when averaged over the volume cor- 
responding to the Beppo Sax field of view. Deeper Hard-X ray 
observations would be required to better compare the two es- 
timates. The values computed here indicate however that they 
can be reconciled. Recently, new hard X-ray observations of 
the Coma cluster have been performed with the new genera- 



A. Bonafede et al.: The Coma cluster magnetic field from Faraday rotation measures 



19 




PQ 



Analytic profile 

Power spectrum 
fluctuations 



i — i *Se-3 : 

A 4 ■■■M. 



2 - 




500 



1000 
r [kpc] 



1500 



2000 



Fig. 20. Left: x 2 plane obtained by comparing simulated and observed cr RM . Right: Profile of the best magnetic field model, magenta 
line refers to the analytic profile (Eq. |6), while the blue line refers to a slice extracted from the simulated magnetic field numerical 
model. Power spectrum fluctuations on the profile are shown. 



tion of satellites (see the work by Wik et al. 2009 using XMM 
and Suzaku data, Lutovinov et al. 2008 using ROSAT, RXTE 
and INTEGRAL data, Ajello et al. 2009 using XMM-Newton, 
Swift/XRT, Chandra and BAT data). These observations failed 
to find statistically significant evidence for non-thermal emission 
in the hard X-ray spectrum of the ICM, which is better described 
by a single or multi-temperature model. Given the large angular 
size of the Coma cluster, if the non-thermal hard X-ray emission 
is more spatially extended than the observed radio halo, both 
Suzaku HXD-PIN and BAT Swift may miss some fraction of the 
emission. These efforts have thus derived lower limits for the 
magnetic field strength, over areas smaller than the radio halo. 
The lower limit reported by Wik et al. (2009) is e.g. (B) > 0.2 
/iG, that is compatible with our results. 



9. Conclusions 

We have presented new VLA observations of seven sources in 
the Coma cluster field at multiple frequencies in the range 1.365 
- 8.465 GHz. The high resolution of these observations has al- 
lowed us to obtained detailed RM images with 0.7 kpc reso- 
lution. The sources were chosen in order to sample different 
lines-of-sight in the Coma cluster in order to constrain the mag- 
netic field profile. We used the numerical approach proposed by 
Murgia et al. (2004) to realize 3D magnetic field models with 
different central intensities and radial slopes, and derived sev- 
eral realizations of the same magnetic field model in order to 
account for any possible effect deriving from the random nature 
of the magnetic field. Simulated RM images were obtained, and 
observational biases such as noise, beam convolution and limited 
sampled regions were all considered in comparing models with 
the data. We found that ctrm and (RM) decrease with increasing 
distance from the cluster center, except for the source 5C4.74, 
that shows a high value of (RM). We argue that this may arise 
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Fig. 21. cr RM and (RM) for the best model (cyan continuous line) 
and its dispersion (cyan dotted lines), given by the rms of the 
different random realizations. Observed points are shown in red. 



from its peculiar position southwest of the Coma cluster core, 
toward the NGC4839 group that is currently merging with the 
Coma cluster. 

Our results can be summarized as follows: 

- the RM ratio and the DP ratio were used to analyze the mag- 
netic field power spectrum. Once a Kolmogorov index is as- 
sumed, the structure-function, the auto-correlation function 
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and the multi-scale statistic of the RM images are best re- 
produced by a model with A max = 34 kpc and A, m „ = 2 kpc. 
We performed a further check to investigate the best value of 
A m ;„ by fitting the Burn law (Burn 1966). This confirmed the 
result obtained from the previous analysis. 

- The magnetic field radial profile was investigated through 
a series of 3D simulations. By comparing the observed and 
simulated cr RM values we find that the best models are in 
the range (B =3.9 nG;i] =0.4) and (B =5.4 fiG;T] =0.7). It 
is interesting to note that the values rj =0.5 and 0.67 lie in 
this range. They correspond to models where the magnetic 
field energy density scales as the gas energy density, or the 
magnetic field is frozen into the gas, respectively. This is ex- 
pected from a theoretical point-of-view since the energy in 
the magnetic component of the intracluster medium is a tiny 
fraction of the thermal energy. Values of Bo >7 fiG and <3 
/iG as well as 77 < 0.2 and 77 > 1 .0 are incompatible with RM 
data at the 99% confidence level. 

- The average magnetic field intensity over a volume of ~ 1 
Mpc 3 is ~ 2 yuG, and can be compared with the equiparti- 
tion estimate derived from the radio halo emission. Although 
based on different assumptions, and although the many 
uncertainties relying under the equipartition estimate, the 
model derived from RM analysis gives an average estimate 
that is compatible with the equipartition estimate. A direct 
comparison with the magnetic field estimate derived from 
the IC emission is more difficult, since the Hard-X detection 
is debated, and depending on the particle energy spectrum, 
the region over which the IC emission arises may change. 
The model derived from RM analysis gives a magnetic field 
estimate that is consistent with the present lower limits ob- 
tained from hard X-ray observations. The values we obtain 
for our best models are still a bit higher when compared with 
the estimate given by Fusco Femiano et al. (2004). It is worth 
to remind, as noted by several authors (see Sec. [8j, that the 
IC estimate derived from Hard X-ray observations could be 
dominated by the outer part of the cluster volume, where 
the magnetic field intensity is lower, depending on the spa- 
tial and energy distribution of the emitting particles. Future 
Hard-X ray missions could help in clarifying this issue. 
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Appendix A: Structure Function and Multi-Scale Statistics with different power spectrum models 

In this Appendix we discuss how other power-law spectral models could be representative of the data presented in the paper. The 
analysis is performed on the basis of the the structure-function, auto-correlation function and multi-scale statistics(MSS). Following 
the approach discussed in Sec. 17.11 we have obtained simulated RM images from different power spectrum models and compared 
them with observed data. We show in Fig. IA. li the structure function, auto-correlation function and MSS for a Kolmogorov power 
spectra that has different A max . We show in Fig. lA.2l fhe fit to the structure, auto-correlation functions and MSS for power spectrum 
models with n = 2, and different values of A max . We show only the plots obtained for the central source 5C4.85. These figures 
demonstrate how the RM data presented in this paper are sensitive to different power spectrum models. 

We note that Kolmogorov power spectra with A max ~100 and 10 kpc fail in reproducing the (RM). These trends can be easily 
understood since power spectrum models with n > 3 have most of the magnetic energy on large spatial scales, and thus small 
changes in A max have a consistent impact on the resulting statistics. According to results presented in Sec. 17.1.11 the case A max =20 
kpc gives a reasonable fit to our data. The best fit is however achieved for A max =34 kpc. In Fig. IA.2l similar fits obtained for power 
spectra models with n = 2 are shown. As indicated by the analysis performed in Sec. 17.1.11 in this case the best agreement with 
observations is achieved for A max of order of hundreds kpc. We note that because of the power spectrum degeneracy, it is possible 
to obtain a reasonable fit to our data. Indeed the case A max =400-800 kpc can reproduce the MSS statistics, although they fail in 
reproducing the S (r) trend on large spatial scales, indicating that a larger value of n is required. 
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Fig. A.l. Fit to the RM images for different Kolmogorov power spectra for the central sources 5C4.85. The different models are 
indicated by different colors (see labels) left: fit to the cr RM and (RM); right: fit to the S(r) and A(r). Red diamonds represent the 
observed statistics. Lines represent the value obtained by averaging 10 power spectra generated with different random seeds. 



Appendix B: Limits on the magnetic field profile from background radio sources. 

Although several arguments (see Sec. 14.2b suggest that the main contribution to the observed RMs is due to the ICM, the best way 
to firmly avoid any kind of local contribution would be to consider only background radio galaxies in the analysis. This is however 
not trivial in general and not feasible here. In fact, sources located in the inner region of the cluster, at distances < (1 — 2)rc are 
fundamental to constrain the magnetic field strength and radial decline, so that a peculiar cluster where background bright and wide 
sources are seen in projection very close to the cluster center would be required. These conditions are not fulfilled in the case of the 
Coma cluster, even though it is a very nearby cluster, where several lines of sight can be inspected. We show in Fig. lB.lK left panel) 
the^ 2 plane obtained by considering only the background radio galaxies: 5C4.74, 5C4.114, 5C4.127, 5C4.42 and 5C4.152. This 
plot shows that for every value of Bo it is possible to find a value of 77 that can reproduce the observed data within l-cr confidence 
level. In the same Fig. in the right panel we show the trends of cr RM and (RM) as a function of r obtained for different value of Bo- 
They all lie within lcr confidence level of the^ 2 plane. It is clear from this plot that the missing information at projected distances 
r < 300 kpc does not permit us to infer the magnetic field strength and radial decline in the Coma cluster. Even unrealistic models 
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Fig.A.2. Fit to the RM images for different power spectra with n — 2 for the central sources 5C4.85. The different models are 
indicated by different colors (see labels) left: fit to the cr RM and (RM); right: fit to the S{r) and A(r). Red diamonds represent the 
observed statistics. Lines represent the value obtained by averaging 10 power spectra generated with different random seeds. 



where r\ <0 cannot be ruled out when the two wide central sources are not considered. Future instruments such as SKA are expected 
to detect many more radio sources seen through a single cluster, possibly allowing this kind of analysis to be performed excluding 
cluster members. 




Fig. B.l. Left: x 2 plane obtained by comparing simulated and observed cr RM for background sources. Lines refer to 1,2 and 3-cr 
confidence level. Right: ctrm and {RM) trends for different models that lie within 1-cr confidence level of the^- 2 . 



